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Abstract 

We will present the result of an analytical calculation of the second order contribution to the 
forward-backward asymmetry A^ B and the shape constant a H for heavy flavour production 
in e + e _ -collisions. The calculation has been carried out by assuming that the quark mass 
is equal to zero. This is a reasonable approximation for the exact second order correction 
for charm and bottom quark production at LEP energies but not for top production at 
future linear colliders. Our result for A^ B is a factor 2.6 (charm) and 4.7 (bottom) larger 
than obtained by a numerical calculation performed earlier in the literature. We study the 
effect of the second order corrections on the above parameters including their dependence 
on the renormalization scale. Further we make a comparison between the fixed pole mass 
and the running mass approach. 
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Experiments carried out at electron positron colliders like LEP and LSD have provided us 
with a wealth of information about the constants ]I[ appearing in the standard model of the 
electroweak and strong interactions. One among them is given by the electroweak mixing angle 
defined by 9w which can be very accurately extracted from the forward-backward asymmetry in 
heavy flavour production in particular when the flavour is represented by the bottom quark ||. 
Recently this quantity is obtained for the charm quark || and in the future one also hopes to 
measure it for the top quark at the large linear e + e~-collider (see e.g. The forward-backward 
asymmetry is extracted from the differential cross section given by 



da H (Q 2 ] 
dcos9 



~(1 + cos 2 9) o vv {Q 2 )r T {p) + a AA (Q 2 )/£(p) 
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(1) 



where 9 is the angle between the outgoing quark H and the incoming electron. The CM energy is 
denoted by Q and the scaling variable p is defined by p = Am 2 /Q 2 where m stands for the quark 
mass. Notice that the form of the cross section above is correct if only final state corrections are 
present which is the case for QCD investigated in this paper. For electroweak corrections (see 
HI) 0) which occur in the initial as well as in the final state, including interference terms, the 
above formula has to be modified. The Born cross sections for quark final states appearing in 
Eq. ([J]) can be written as 
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where N denotes the number of colours in the case of the gauge group SU(N) (in QCD one 
has N — 3). Furthermore in the expressions above we adopt for the Z-propagator the energy 
independent width approximation 



Z(Q 2 ) =Q 2 -M 2 z + iM z T z . 



(5) 



The charges of the lepton and the quark are given by and e q respectively and the angle 9y/ 
defined at the beginning appears in the electroweak constants as follows. 
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Cv,u = Ca/ (1 — — sin 2 9w), Cy,d 
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The functions fh (k = T,L,A; I = v, a) in Eq. (IT]) can be computed order by order in perturbative 

LNS 



QCD from the non-singlet quark coefficient function C k 
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where p stands for the factorization as well as the renormalization scale. The quark coefficient 
function appears in the fragmentation function Fk(x,Q 2 ) with x = 2p.q/Q 2 where p is the mo- 
mentum of the outgoing hadron which originates from the quark. These fragmentation functions 
describe the production of the quark and its subsequent decay into a hadron. The forward back- 
ward asymmetry, denoted by Ap B , appears when we divide the expression in Eq. (UJ) by the total 
cross section. The ratio can then be expressed in the following way 
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Further the shape coefficient a H is defined by 
a H (Q 2 ) 
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and the total cross section for heavy flavour production is equal to 

C(Q 2 ) = ^VV(Q 2 ) [frip) + fl(p)] + CTAA(Q 2 ) [ft{p) + flip) • 

The functions f l k can be expanded in the strong coupling constant a s as follows 
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The lowest order contributions corresponding to the Born reaction 

V ^H + H, 

with V — 7, Z are given by 
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where p is defined below Eq. (|IJ). The next-to- leading order (NLO) contributions originate from 
the one-loop virtual corrections to the Born reaction (|13|) and the gluon bremsstrahlungs process 



V ^E + E + g . 



(15) 
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The NLO contributions have been calculated by several groups in the literature for the case 
m 7^ (see [T7I— | TT5|) . However there were some discrepancies between the results so that we 
decided to calculate them again using a different method. Here they are derived from the order 
a s contribution to the coefficient functions for heavy quarks which appears in the integrand of 
Eq. (0). Our computations lead to the same answers as given in appendix A of [11]]. After 
performing the integral in Eq. (0) we obtain the following results for k = T,L 
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(20) 



l + v ^7 

Further the colour factor Cp is given by Cp = (N 2 — 1)/2N. The functions appearing 
above are defined by 
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(23) 
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where ({n), which appears for n = 2, 3 in the formulae of this paper, represents the Riemann 
^-function and Li 2 (x) denotes the dilogarithm. Using Eqs. dl6|) -(|l9|) one can check that the order 
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a s contribution to the total cross section in (|TT ) is in agreement with the literature |12| (for the 
vector part see also [13]). Integration of the asymmetry quark coefficient function provides us 
with the result 
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The functions (k — T,L; I = v,a) are related to the functions if 2 , and ii/g presented in 
Eq. (15) and appendix A of |§ and we agree with their result. The same holds for which 
is proportional to H5 in the reference above. There is also agreement with the calculation in 
||10|| . The comparison is made by expanding the functions above and those in || up to seven 
powers in p. The next-to-next-to-leading order (NNLO) contributions come from the following 
processes. First one has to compute the two-loop vertex corrections to the Born process ( |13"D 
and the one-loop corrections to (pip. Second one has to add the radiative corrections due to the 
following reactions 



V -^E + E + g + g. 
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H + H + H + H. 
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V -> E + E + q + q, 



where q denotes the light quarks. The results for f 1 ^ presented below are computed for the 
contributions where the vector boson V is always coupled to the heavy quark H so that the 



(30) 



light quarks in Eq. (pCf ) are only produced via fermion pair production emerging from gluon 
splitting. Besides these contributions there are other ones which have been treated in ||14|| . The 
latter consist of all one- and two-loop vertex corrections which contain the triangular quark-loop 
graphs. Following the notation in Jill then contribution to the forward-backward symmetry 
will be denoted by Fq^^ and Fq~qd respectively. They only show up if the quarks are massive 
and are coupled to the Z-boson via the axial-vector vertex. Notice that one has to sum over all 
members of one quark family in order to cancel the anomaly. Adopting the mass assignment in 
[HI we take the top to be massive and put the other quark masses, including that of the bottom, 
equal to zero. Further in |TJ| one has included all terms originating from reaction (|30| ) where 
diagrams with light quarks attached to the vector boson V interfere with those describing the 
coupling of the heavy quarks to the vector boson. Notice that this contribution denoted in 
by Fq CD vanishes if all quarks including H are taken to be massless provided one sums over all 
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members in one family. However we will omit that part of reaction ([30]) where the heavy quarks 
are produced via gluon splitting and the light quarks q are coupled to the vector boson V. This 
contribution denoted by Fq^ co in [13] needs a cut on the invariant mass of the heavy flavour pair 
and it was computed for the first time in JT5|]). Finally notice that these additional contributions, 
denoted by Fqcd above, only show up in order a 2 . Moreover if we put m = they only appear 
in the forward-backward asymmetry @ but cancel between numerator and denominator in the 

/ (2) 

shape coefficient flTU|). The results for / fe follow from the transverse and longitudinal coefficient 
functions calculated in |16[ and the asymmetry coefficient function computed in |17| . Because of 
the complexity of the calculation of these functions the heavy quark mass was taken to be zero. 
This approximation is good for the charm and bottom quark but not for the top quark as we 
will see below. Substituting these coefficient functions in the integrand of Eq. (0) we obtain 
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Here the colour factors are given by Ca = N and Tf = 1/2 (for Cp see below Eq. 
Further rif denotes the number of light flavours which originate from process (PUp. Finally \i 
appearing in the strong coupling constant a s and the logarithms in Eqs.(|3T|), (p^) represents 
the renormalization scale. Notice that the coefficient of the logarithm is proportional to the 
lowest order coefficient of the /3-function. The same holds for ((3) in Eqs. ([31]) and d33|). The 



logarithm does not appear in f A 
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zero there is no distinction anymore between f^'^ and f^'^ (k 



Since the mass is equal to 
T, L) unlike in the case for 
the first order corrections in Eqs. (|TED-(|ll?D where the heavy quark was taken to be massive. 
Furthermore one can check that substitution of Eqs. (|31|), (|32|) into ([□]) provides us with the 
order a 2 s contribution to the total cross section which is in agreement with the results obtained 
18]. For zero mass quarks the forward-backward asymmetry becomes equal to 
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For the discussion below and the notation often used in the literature (see 
to write Ap B in the following way 




^b(Q 2 )=^pb (0) (^ 



0L 8 (H) 
1 Ci 



71 



71 



C 2 



(35) 



5 



The order a 2 contribution presented above has been also calculated in [14| but then in a numerical 
way. Unfortunately we get a different answer. First we disagree with the statement above Eq. 
(48) in [TJj that reaction (|29"D does not contribute to Ap B . For the latter we obtain the following 
contribution to C2, denoted by A^C2, which equals 



A«c 2 = -(Cf - \CaCf)^ 
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0.14 



(36) 



The result for reactions ( F2S| ) and (|3T)D, including the virtual corrections, can be obtained by 
subtracting Eq. (|3^) from Eq. fl3"4"|). Using the notation in Eq. (45) |14| and choosing nf = 5 
we find the following contributions 
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-2.83 iViV c = -42.38 TT R = 13.75 . 
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From the results above we infer that the discrepancies mainly occur in the Cp and C^CV-terms 
of Eq. (|34"D . Substituting Cp = 4/3 in the expression above we obtain for the coefficient of the 
(a s (Q)/7r) 2 -term the value —9.49 instead of —2.6 quoted in Eq. (4.7) of |TJ| which amounts to 
a discrepancy of a factor of about 3.7. Notice that the bulk of the second order correction to 
Eq ( P4"|) is coming from f% in Eq. (1331) which amounts to 9.216. The remaining part can be 
traced back to the functions / fe in Eqs. ([H]), (0) leading to the contribution 0.409. 



Finally we are now also able to present the second order correction for the shape coefficient (10) 
in an analytical way 
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which can be written in the same way as has been done for A^ B in Eq. (BH) so that we get 
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We will now discuss the effect of the order ot 2 s corrections on the forward backward symmetry and 
the shape coefficient for the bottom and the charm quark. Further we omit any effect coming 
from the electroweak sector. Our results are obtained by choosing the following parameters 
(see pi). The electroweak constants are: M z = 91.187 GeV/c 2 , T z = 2.490 GeV/c 2 and 
sin 2 9w = 0.23116. For the strong parameters we choose : A^ IS = 237 MeV/c (jif = 5) which 
implies a s (Mz) = 0.119 (two-loop corrected running coupling constant). Further we take for 
the renormalization scale fi = Q unless mentioned otherwise. Notice that we study Ap B and a H 
for H = c, b at the CM energy Q = Mz- For the heavy flavour masses the following values are 
adopted : m c = 1.50 GeV/c 2 , m b = 4.50 GeV/c 2 and m t = 173.8 GeV/c 2 . The results for the 
bottom quark can be found in table |l[ The values for c 2 are obtained by adding to Eq. (|34D the 
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^FB 




m b = 4.50 GeV/c 2 


m b (M z ) = 2.80 GeV/c 2 


Cl 


0.789 


0.858 


c 2 


8.89 (9.63) 


8.89 (9.63) 


a(o) 

^FB 


0.1052 


0.1052 


^FB 


0.1020 (- 3.04 %) 


0.1017 (- 3.33 %) 


^FB 


0.1007 (- 4.28 %) 


0.1004 (- 4.56 %) 




di 


2.08 


2.33 


d 2 


23.0 


23.0 


a<°> 


0.994 


0.998 


a« 


0.916 (- 7.85 %) 


0.910 (- 8.82 %) 


a( 2 ) 


0.883 (- 11.2 %) 


0.877 (- 12.1 %) 



Table 1: The forward-backward asymmetry and the shape constant of the bottom quark. 



contributions mentioned below Eq. (|30|). They were denoted by Fqc^ ', Fqgd an d Fqcd 1X1 EI 
For the value of the top mass given above we obtain for bottom production c^qcd = —0.645, 
c 2~qcd = —0.218 and c^qcd = 0.123. In the entry for c 2 we have also mentioned between the 



brackets the result obtained from the contribution of (|3~4] ) which is equal to c 2 = 77/8. Notice 
that the second order contributions to the quantities in the tables are obtained by multiplying 
Apg at m = with c 2 . Furthermore we have put in the tables the correction in percentages 
of the radiatively corrected quantities Ap B , a H with respect to their zeroth order result. From 
table [l] we infer that the order a s as well as the order a 2 s corrections to both Ap B and a b are 
negative. In the case of the forward-backward asymmetry the QCD corrections are moderate in 
particular the second order ones. The latter would be even smaller if we had taken the value 
c 2 = 1.9 quoted in |TJ] which is a factor 4.7 less with respect to our result. In the case of 
the shape coefficient the QCD corrections are at least twice as large. The latter are reduced 
if for the reference axis the thrust axis is taken instead of the quark axis 0. Notice that the 
quark axis has been chosen in our calculation. Further we want to mention that the zero mass 
approximation for the second order coefficient c 2 is quite reasonable. This is revealed by the first 
order coefficient when the quark mass is chosen to be zero which leads to the values c\ = 1 Eq. 
(|34]) and d\ = 8/3 Eq. fl38D . In the case of the bottom quark we observe a deviation of 11 % for 
Ci whereas for d\ it amounts to 22 % which is twice as large. For the charm quark (see table 0) 
these values become smaller i.e. about 7.5 % for both coefficients. If we expect that the same 
deviations occur for the coefficients c 2 Eq. (|34])and d 2 Eq. ( |38"D one gets a reasonable estimate of 
the theoretical error on the second order corrections. We also studied the effect of the running 
quark mass on the forward-backward asymmetry and the shape coefficient. For this purpose one 
has to change the on-mass shell scheme used in Eqs. (|i~6D - (|i~9"D , (p5|) into the MS-scheme. This 
can be done by substituting in all expressions the fixed pole mass m by the running mass . 
Moreover one has to add to the first order contributions (|16|)-(fL9), (|25|) the finite counter term 



4 - 3 In 



m 2 (/i) 



'dff\ P ) 



d m 



(40) 

m=m(fi) 



7 



^FB 




m c = 1.50 GeV/c 2 


m c (M z ) = 0.662 GeV/c 2 


Cl 


0.924 


0.962 


c 2 


11.5 (9.63) 


11.5 (9.63) 


AO) 

^FB 


0.0750 


0.0750 


^FB 


0.0724 (- 3.47 %) 


0.0723 (- 3.60 %) 


4(2) 
^FB 


0.0712 (- 5.07 %) 


0.0710 (- 5.20 %) 


a c 


d\ 


2.46 


2.58 


d 2 


23.0 


23.0 


a(°) 


1.000 


1.000 


a« 


0.907 (- 9.30 %) 


0.903 (- 9.70 %) 


a( 2 ) 


0.874 (- 12.6 %) 


0.870 (- 13.0 %) 



Table 2: The forward-backward asymmetry and the shape constant of the charm quark. 



where /i stands for the mass renormalization scale for which we choose \x = Q. Further we adopt 
the two- loop corrected running mass with the initial condition m(/i ) = /^o- Using the relation 
between the MS'-mass and the fixed pole mass, as is indicated by the first factor on the right- 
hand side in Eq. fl4"0|), we have taken for bottom production /xq = 4.10 GeV/c 2 which corresponds 
with a pole mass nib = 4.50 GeV/c 2 . This choice leads to mj(M^) = 2.80 GeV/c 2 which is 5 % 
above the experimental value 2.67 GeV/c 2 measured at LEP pp| . The results are presented in 
the second column of table |l|. Comparing the latter with the first column we observe that the 
QCD corrections become a little bit more negative. The values of Apg (i = 0, 1,2) decrease a 
little too. The same features are also shown by except for which slightly increases. 

In table |2| we also present results for the charm quark. In the case of the running mass we 
have chosen /i = 1.30 GeV/c 2 so that the pole mass becomes m c = 1.50 GeV/c 2 . This leads to a 
value of m c (Mz) = 0.662 GeV/c 2 which is rather low. Furthermore for charm quark production 
the additional contributions become c^qcd = 1-359, (^qcd = 0-323 and cf q CD = 0.211. Our 



result for c 2 is about 2.6 times larger than the value 4.4 obtained in [ 14"| . The features are 
the same as for the bottom quark. The differences between the numbers in the left-hand and 
right-hand column in table |^ become even less. The reason that the running of the mass in the 
case of the charm and the bottom quark hardly introduces any effect on the forward-backward 
asymmetry and the shape constant can be attributed to the fact that the mass of both quarks 
are small with respect to the CM energy so that the phase space is not much affected. Moreover, 
as far as the dynamics is concerned, these quantities are not proportional to the mass. This is 
also revealed by the constants c\ and d\ which do not deviate very much from their zero mass 
values. These arguments do not apply to the top quark except when Q 3> m t . This is shown 
in table ^ where we have studied the above quantities at a CM energy Q = 500 GeV/c. For 
the running mass we have chosen /i = 166.1 GeV/c 2 so that the pole mass becomes equal to 
m t = 173.8 GeV/c 2 . This leads to a value m t (Q) = 153.5 GeV/c 2 . The constants c\ and d\ in 
the perturbation series completely differ from the ones given at m = at which they become 
1 and 8/3 respectively. Therefore the running mass will have a large effects on these constants 
which is revealed by the switch of sign in table Hence the Born approximations to Ap B and 
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A 1 

^FB 




m t = 173.8 GeV/c 2 


m t (Q) = 153.5 GeV/c 2 


Cl 


- 0.757 


2.11 


4(0) 
^FB 


0.407 


0.455 


^FB 


0.417 (2.46 %) 


0.426 (- 6.37 %) 


a 1 




- 0.911 


5.13 




0.406 


0.514 




0.417 ( 2.71 %) 


0.435 (- 15.4 %) 



Table 3: The forward-backward asymmetry and the shape constant of the top quark at Q — 
500 GeV/c. 

a* will change while going from the fixed pole mass to the running mass approach. However the 
order a s corrected quantities are less sensitive to the choice between the running or the fixed 
pole mass because of the compensating term in Eq. fl40|) . From the above it is clear that the zero 
mass approximation to C2 and 62 makes no sense in case of the top quark and we have omitted 
these contributions to Ap B and a* in table || Finally we want to comment on the renormalization 
scale dependencies of the forward backward asymmetry and the shape constant. If we vary the 
scale /i between Q/2 and 2 Q the changes in Ap B are small. It introduces an error of 0.002 for 
the bottom quark and 0.003 for the charm quark. For one can draw the same conclusion and 
the error becomes 0.005 and 0.007 respectively. In the case of the top quark a variation in the 
renormalization scale makes no sense because of the missing order a 2 s correction. Its computation 
for massive quarks will be a enormous enterprise. 

Summarizing our findings we have computed the order a 2 contributions to the forward- 
backward asymmetry and the shape constant in an analytical way provided the heavy flavour 
mass is chosen to be zero. Further we found a discrepancy with a numerical result calculated 
earlier in the literature for Api? 2 ''- ^he secon d order corrections are noticeable. The transition 
from the fixed pole mass to the running mass approach does not introduce large changes in the 
values of v4p B and a H except for the first order constant in the perturbation series when H = t. 
This indicates that the zero mass approach breaks down unless Q 3> m. Also a variation of the 
renormalization scale does not lead to large effects. The latter are almost equal to the differences 
between the results obtained by the fixed pole mass and the running mass approach. 
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